Folia Forestalia PolonUa, series A, 2014, Vol. 56 (2), 93-104 



ORIGINAL ARTICLE 



DOI: 10.2478/ffp-2014-0009 



Height-diameter models for mixed-species forests 
consisting of spruce, fir, and beech 

Rudolf Petras 1 El, Michal Bosel’a 2 , Julian Mecko 2 , Julius Oszlanyi 3 , Ionel Popa 4 

1 National Forest Centre, Forest Research Institute Zvolen, T.G. Masaryka 22, 960 92 Zvolen, Slovakia, 
phone: +421 455314231, e-mail: petras@nlcsk.org 

2 National Forest Centre, Forest Research Institute Zvolen, T.G. Masaryka 22, 960 92 Zvolen, Slovakia 

3 Slovak Academy of Sciences, Institute of Landscape Ecology, Stefanikova 3, 814 99 Bratislava, Slovakia 

4 Forest Research and Management Institute, Forest Research Station for Norway Spruce Silviculture, 
Calea Bucovinei 73 bis, 725100 Campulung Moldovenesc, Romania 



Abstract 

Height-diameter models define the general relationship between the tree height and diameter at each growth stage 
of the forest stand. This paper presents generalized height-diameter models for mixed-species forest stands consist- 
ing of Norway spruce (Picea abies Karst.), Silver fir ( Abies alba L.), and European beech ( Fagus sylvatica L.) from 
Slovakia. The models were derived using two growth functions from the exponential family: the two-parameter Mi- 
chailoff and three-parameter Korf functions. Generalized height-diameter functions must normally be constrained 
to pass through the mean stand diameter and height, and then the final growth model has only one or two parameters 
to be estimated. These “free” parameters are then expressed over the quadratic mean diameter, height and stand age 
and the final mathematical form of the model is obtained. The study material included 50 long-term experimental 
plots located in the Western Carpathians. The plots were established 40-50 years ago and have been repeatedly 
measured at 5 to 10-year intervals. The dataset includes 7,950 height measurements of spruce, 21,661 of fir and 5,794 
of beech. As many as 9 regression models were derived for each species. Although the “goodness of fit” of all models 
showed that they were generally well suited for the data, the best results were obtained for silver fir. The coefficient 
of determination ranged from 0.946 to 0.948, RMSE (m) was in the interval 1.94-1.97 and the bias (m) was -0.031 to 
0.063. Although slightly imprecise parameter estimation was established for spruce, the estimations of the regression 
parameters obtained for beech were quite less precise. The coefficient of determination for beech was 0.854-0.860, 
RMSE (m) 2.67-2.72, and the bias (m) ranged from -0.144 to -0.056. The majority of models using Korf’s formula 
produced slightly better estimations than Michailoff’s, and it proved immaterial which estimated parameter was 
fixed and which parameters were free. 
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Introduction 

Height growth of trees is their one of the most charac- 
teristic biological feature. It is generally accepted that 
tree height depends not only on age, diameter, species, 
and species mixture but also on the quality of the site 
where the trees grow. Height is usually defined by spe- 
cific models, where the tree height h is the function of 
its diameter d\ 

h = f(d) (1) 

While individual height-diameter models were ini- 
tially defined using simple mathematical functions, such 
as polynomial function (Naslund 1929; Assmann 1943; 
Petterson 1955; Kennel 1965), fractional polynomials 
(Korsun 1935; Michailoff 1943) or into the exponential 
functions (Korf 1939; Freese 1964; Sharma and Parton 
2007). More complex functions had to be transformed 
into linear form (Naslund 1936; Prodan 1944). Schmidt 
(1967) tested 6 functions and the best results were ob- 
tained by Korsun’s function. Kennel (1972) tested the 
same functions in beech stands in Bavaria and Swit- 
zerland. Meanwhile, Sterba and Marschall (1976) tested 
8 functions using data from the National Forest Inven- 
tory in Austria, and they reported that several functions 
gave quite satisfactory estimations. Zeide (1989, 1993) 
analyzed many popular growth equations and claimed 
that Korf’s (1939) equation is substantially more accu- 
rate than other growth equations. Its standard error of 
estimate was 2. 1-4. 8 times less than the errors of the 
Chapman-Richrds, Weibul, Gompertz, and logistic 
equations, respectively. 

Height-diameter curves have a special shape and 
position in each growth stage of forest stands, and these 
emanate from stage-dependent height-diameter mod- 
els. In this sense, the steeper curve of the older stand is 
placed above that of the younger one. It has an elongated 
interval towards larger-diameter trees, but this is short- 
er for small-diameter trees. However, the height-diame- 
ter curve-shift between growth stages is quite irregular 
during stand development. It is lower for small-diame- 
ter trees in older stands and higher for large-diameter 
trees in younger stands. These properties of height-di- 
ameter models can be generalized in one model, where 
the tree height h is a function of its diameter d and other 
stand variables such as either the mean diameter d m and 



height h m , together with the stand age t, are defined in 
the following equation: 

h = f(d,d m ,h m ,t ) (2) 

The models using Equation (2) were constructed 
sequentially. The models were first fitted by the known 
functions in formula (1) and these parameters were sub- 
sequently estimated using selected stand characteristics 
as independent variables. One parameter was omitted 
from fitting and later employed in adjusting the final 
mathematical formula so that the curves passed through 
the mean stand diameter and height (d m , h 11L ). Kennel 
(1972) used equation (2) to derive a model from 381 
height-diameter curves of beech stands for the 3 thin- 
ning intensities - light, moderate, and intensive. For this 
purpose, he utilized the function provided by Petterson 
(1955). Smelko et al. (1987) derived height-diameter 
models for 12 species using the function published by 
Michailoff (1943). In other countries, the function is 
known as Schumacher's (1939), and it was developed on 
the basis of graphical models published by Halaj (1955). 
Nagel (1991) derived two models for red oak. One was 
based on Kennel (1972) and the second on the Sloboda et 
al. (1993) function, based on Michailoff (1943). Similar to 
Smelko et al. (1987), he fitted the regression parameters 
using mean stand diameter and height as independent 
variables. Using this same formula, Petras and Mecko 
(2005) derived their height- diameter models for poplar 
clones. After testing 3 height-diameter models, Hui and 
Gadow (1993) proposed their own model which consist- 
ed of three allometric functions. Therein, tree heights 
are modelled as a function of their diameter and top 
stand height. In addition, Temesgen and Gadow (2004) 
developed a generalized height-diameter model using 
a database of permanent sample plots in British Colum- 
bia. This is based on a Weibull-type function, where the 
relationship between diameter and height is influenced 
by the relative competitive position of the trees and by 
stand-density variables. Stankova and Dieguez-Aranda 
(2013) tested ten generalized and six local height-diame- 
ter models for pine forests in Bulgaria. 

The tree heights in forest stands serve as the ba- 
sis for evaluation of their growth and production. This 
is especially essential in mixed-species forests such as 
spruce-fir-beech, where these three different growth- 
type species grow in mutual competition. The objec- 
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tives of this paper is to develop (i) local height-diameter 
models and forest stand characteristics for the species; 
(ii) global height-growth models for spruce-fir-beech 
mixed forests; and (iii) to evaluate height variability, 
precision and accuracy of the models using a long term 
experimental data for spruce-fir-beech forests. 



Material and methods 

Study site and empirical data 

Empirical data used in this study included 252 repeat- 
ed measurements in 50 long-term experimental plots 
(LTP) established in the 1960’s and 1970’s in the West- 
ern Carpathians in Slovakia. The plots were established 
to study the growth and production of homogeneous 
and mixed-species forest stands, and they are located 
in the Pol’ana Massif (48°36'43"N and 19°36'37''E) 
and the Slovenske Rudohorie Mts. (48°50'20"N and 
20°44'06"E) at altitudes ranging from 480 to 970 m. 
The sites were classified in Abieto-Fagetum, Fageto- 
Abietum, Fagetum abietino-piceosum, and Fagetum 
quercino-abietinum forest types, after Zlatnik (1959, 
1976). The species mixture varied in the plots. All 
three investigated species were present in 22 LTP, while 
23 LTP had spruce and fir admixture, 2 LTP had spruce 
and beech, and the remaining 3 LTP had fir and beech. 
Fir was most prominent, followed by spruce and then 
beech. All the study forest stands were one-layered and 
even-aged. The stand-age at the plots ranged from 32 
to 159 years, and stands were repeatedly measured and 
treated by light thinning; mainly at 5-year intervals. 
The majority of the plots were measured 4-8 times, 
and the stand age ranged from 73 to 202 years after 
a 40-year period. The area of LTP ranged from 0.200 
to 1.215 ha. LTP size, when established, was defined 
so that it included at least 300 trees. All trees in the 
plots were numbered, marked and had their diameter 
at breast-height (DBH) and absolute height measured. 
Diameter was measured using a calliper with precision 
of 1 mm for the whole study period, however, height 
was measured by different hypsometers. First it was 
Blume-Leiss, later Spiegel-Relaskop and since 2002 it 
was ultrasonic VERTEX. In addition, the trees’ social 
position and both the stem quality and damage were as- 
sessed. Although the height of all trees was measured 
only on the first and last experimental occasions, tree- 



height was also measured between these periods for the 
trees sampled for height-diameter curve construction. 

Model construction 

The height curve from equation (1) formed the basis for 
the general height-diameter model in equation (2). After 
thorough review of published results (Zeide 1989, 1993; 
Mehtatalo 2004) and the author’s own experience, two 
functions were selected for the model. 

Firstly Michailoff’s (1943) function formed equa- 
tion (3): 

f-1 

h(d) = 1.3+ a.e ' d (3) 



where d denotes the diameter at breast height, and a and 
b are the regression parameters to be estimated. 



This function has several properties. It has an as- 
ymptote, with inflexion point and high flexibility. It 
is quite simple and has the following two relatively 
stable parameters; a is the intercept which the largest 
tree heights approximate and b defines the shape of the 
height-diameter curve. 

Secondly, Korf’s (1939) function formed equa- 
tion (4): 



h(d) = 1.3+ A - exp 



I — n 



-■d [ 



■ -») 



(4) 



This function has three parameters (A, k, and n) 
and, according to the author’s experience, it is much 
more flexible than the previous one. 

To achieve the requirement that each model curve 
must pass through the mean stand diameter and height 
(d m , h,J (Smelko et al. 1987; Nagel 1991; Sloboda et 
al. 1993; Petras and Mecko 2005), one parameter must 
be fixed while the others remain free and are later es- 
timated depending on particular stand variables. Two 
variants were selected for equation (3). The first had pa- 
rameter a fixed and b free, depending on the mean stand 
diameter and height (d q , h q ). This model has the simpler 
mathematical form (Smelko et al. 1987): 



, , / , \ 






h{d,d q ,hF) = \3 + \h q -1.3J-exp 


b ■ 


d d 

\ i / J 



where d q , h q - mean stand quadratic diameter and 
height. 
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The variability of the parameter b was analyzed 
in relation with mean stand diameter dq , mean stand 
height hq, and stand age t. Product and sum of two pow- 
er functions and an exponential function were found as 
the most appropriate functions for this purpose. Three 
variants were tested for parameter b\ 

- Michailoff 1 

b(d q ,h q ) = pl + P 2-h p3 -d p4 (5.1) 

- Michailoff 2 

b{d q ,h q ) = p\ + p2- h q + p2> ■ d q (5.2) 

- Michailoff 3 

b(t) = p\ + p2- e 1 (5.3) 

where: t is the stand age. 



Similar to the Michailoff’s function, the following 
equations were tested for parameters k and n: 

- Korfl 



k(h ,d ) = pl-h pl -d p 3 

v q q' r q q 

«( h q ,d q ) = p4- h p5 ■ d p6 

- Korf 2 

p2 

k{d q ) = p\-e dq 

p4 

n(d q ) = p2-e dq 

- Korf 3 



k(t) = pi ■ t p2 
n(t) = p2-1 pl[ 



(7.1) 



(7.2) 



(7.3) 



where parameter k is fixed, and A and n are free, the 
model has the form: 



The second variant of Michailoff’s function had pa- 
rameter a free and b fixed; in the form: 



h( d,d q ,h q ) = \3+ a - exp 



-d 



-In 



h -1.3 

V i 7 



( 6 ) 



where In is the natural logarithm. 



Two variants were tested for parameter a: 
- Michailoff 4 



a(d q ,h q ) = pi- h p2 ■ d p3 (6. 1) 

- Michailoff 5 

a{d q , h q ) = pi + p2 ■ h p q ■ d pA (6.2) 



Korf’s function (4) has three parameters; one is 
fixed with an implicit condition while the remaining 
two are free and depend on individual stand variables. 
Where parameter A is fixed and parameters k and n are 
free, the model has the following form: 



h(d,d q ,h q ) = 1.3 + 






•exp 



1 — n 



■(d^-d^-f 



(7) 



h(d,d q ,h q ) = 1.3 + 



+^4-exp 



In 



( , , / y 

*- — 1-3 1 d 



v A 



d 

V 9 7 



( 8 ) 



where In is the natural logarithm. 

The following equations were tested for parameters 
A and n: 

- Korf 4 



A(d ,h ) = pl-h p2 -d p3 

v q q y c q q 

n(d ,h ) = p4 ■ h pi ■ d pb 

v a ' a ' r a a 



( 8 . 1 ) 



- Korf 5 



A(d q ,h q )= pl + p2-hfd p4 
n(d ,h) = P 5-b p6 -d p7 

x a ' a ' r a a 



( 8 . 2 ) 



- Korf 6 



A(d q ,h q ) = pl + p2h p3 -d p4 
l — n(d ,h ) = p5 + p6- h pl ■ d pS 



(8.3) 



The third variant of the Korf function, where the 
parameter n would be fixed and the other parameters 
free was not tested, because the parameter n best de- 
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termines the shape of the curve. It is desired to have 
the parameter n free. The mathematical form of the 
generalized height-diameter models according to for- 
mula (2) were derived by adding the equations for the 
free parameters to the general formula. Here, equations 
(5.1)— (5.3) were added to model (5), (6.1) and (6.2) to 
formula (6), (7.1) — (7.3) to (7), and (8.1)— (8.3) to (8). 



Heights at the long-term experimental plots 

The height-diameter models for each repeated measure- 
ment and species were derived using Michailoff’s func- 
tion (3). Parameters a and b were tested and fitted using 
the stand age t as the independent variable. The following 
allometric functions were selected as the most relevant: 

a(t) = a t ■ t al (9) 

b(J)=b ] -t h2 (10) 



The allometric function is simple and has only 
2 parameters, compared to the one proposed by Johann 
(1990). Addition of the functions to formula (3), where 
tree height h is a continuous function of the diameter d 
and age t, provides a model in the following form: 



h(d,t) = 1.3+0, t° 2 



f 

•exp 

v 



-b, • t 



b2 \ 



(ll) 



Regression models were derived for each species 
and each LTP from formula (11), with parameters a h 
a 2 and b h b 2 estimated using the Gauss-Newton opti- 
mization technique in a non-linear least-squares proce- 
dure using QC. Expert (Kupka 2008). The precision and 
accuracy of the models were then evaluated using root 
mean square error (RMSE), coefficient of determina- 
tion (R 2 ), and bias. 

The quadratic mean stand diameter d q for each 
species and each repeated measurement was calcu- 



lated for each LTP. The mean stand height h q was cal- 
culated using equation (11). Hence, the (5 x n) matrix 
was created for each species from experimental data. 
This comprised 5 columns where mutually related 
tree characteristics (d, h), mean stand variables (dq, 
hq), and stand age t were listed. The row number (n) 
presents the number of trees with measured heights. 
There were 21,661 silver fir, 7,950 spruce and 5,794 
beech height measurements. 



Results 

Height-diameter models at the long-term 
experimental plots 

The height-diameter equations (using equation 11) were 
derived for individual tree species at the LTP. The coef- 
ficient of determination (R 2 ) for all three species (tab. 1) 
ranged from 0.555 to 0.970 and the root mean square 
errors (RMSE) were from 1.07 to 3.72 m. Their arith- 
metic averages had variation; R 2 = 0.815-0.870, and 
RMSE = 1.99 m for fir, 2.02 m for spruce and 2.49 m for 
beech. The highest variability and RMSE were in LTP’s 
where low numbers of trees were measured. 

The reason for such a high error rate was the wide 
variability in tree-height in individual forest stands, 
but some errors may also be due to imprecise height 
measurement. This is exemplified by the height-diam- 
eter models of silver fir from repeated measurements 
on LTP 83 (fig. 1). Here, differences in tree height at 
approximately the same diameter ranged from 3 to 
5 m. The model precision of RMSE was 1.27 m and 
R 2 = 0.962. However, there was obvious bias between 
the model and measurements at stand ages of 46 and 
79 years, and we assume that this was due to imprecise 
height measurement rather than the function accuracy. 
Such phenomena are quite frequent in long-term plots 



Table 1 . “Goodness of fit” of the height-diameter function using modified Michailoff function (equation 11) 



Species 


Number 
of LTP 


Number of heights 
per LTP 


R 2 -Coefficient 
of determination 


RMSE-Root 
mean square error (m) 


min-max 


mean 


min-max 


mean 


min-max 


mean 


Norway spruce 


47 


16-856 


177 


0.555-0.970 


0.826 


1.07-3.72 


2.02 


Silver fir 


48 


34-1,098 


470 


0.701-0.962 


0.870 


1.22-3.42 


1.99 


European beech 


27 


23-1,177 


246 


0.601-0.937 


0.815 


1.45-3.31 


2.49 



Folia Forestalia Polonica, series A, 2014, Vo I. 56 (2), 93-104 



98 



Rudolf Petras, Michal Bosel’a, Julian Mecko, Julius Oszlanyi, Ionel Popa 



(Kennel 1972, Johann 1990). Where imprecise meas- 
urements are less frequent and their identification is 
possible, it is better to exclude them from modelling 
or correct them before modelling. In addition, if the 
model curves are not regularly arranged above one an- 
other and it is uncertain which one is wrong, it is better 
to derive them in one generalized model (Curtis 1967; 
Pretzsch 2009). 




Figure 1. Selected generalized height-diameter models 
for silver fir based on repeated measurements on LTP 83, 
at stand ages of 36 years (the first measurement), 46 and 
79 years (the last measurement). R 2 = 0.962, RMSE = 1.27 m 



Species-specific generalized height-diameter 
models 

Although all models were reliable according to ”good- 
ness-of-fit” statistics (tab. 2), the most precise estima- 
tion was for silver fir. Coefficients of determination 
ranged from 0.946 to 0.948, RMSE (m) 1.94 to 1.97 
and bias (m) -0.031 to 0.063. The estimation of the 
parameters for spruce was slightly less precise with 
an R 2 range of 0.921-0.924, RMSE (m) 2.05-2.12, 
and bias (m) -0.007 to 0.117. The precision and ac- 
curacy of the estimation for beech was as follows: R 2 
from 0.854-0.860, RMSE (m) 2.67-2.72 and bias (m) 
from -0.228 to -0.056, which is essentially less pre- 
cise compared to the previous species. The majority 
of the models using Korf’s function gave slightly bet- 
ter estimates than those from Michailoff’s function, 
irrespective of which parameters were fixed or free. 
Although the Michailoff (3) and Korf (3) models had 
free parameters explained only in the relation to stand 
age, they are inferior to other models because they 
have an extra variable to be estimated. The residuals 
(. h-h mod ) greatly varied for each species (fig. 2-4). For 
the all species the residuals were ±10 m. However, the 
majority was up to ±5 m. The largest residuals were 
found for trees with DBH between 15 and 40 cm, and 



Table 2. Statistical characteristics of the generalized height-diameter models for the three species (gray cell indicates 
the smallest sum of squares) 



Model 


Norway spruce 


Silver fir 


European beech 


sum of 
squares 


R 2 


RMSE 

(m) 


bias (m) 


sum of 
squares 


R 2 


RMSE 

(m) 


bias (m) 


sum of 
squares 


R 2 


RMSE 

(m) 


bias (m) 


Michaj 1 


34,551 


0.922 


2.085 


-0.003 


82,769 


0.947 


1.955 


-0.012 


41,956 


0.857 


2.692 


-0.103 


Michaj 2 


34,956 


0.921 


2.097 


-0.006 


84,382 


0.946 


1.974 


-0.015 


42,043 


0.857 


2.694 


-0.100 


Michaj 3 


34,130 


0.923 


2.072 


-0.007 


81,954 


0.948 


1.945 


-0.031 


41,284 


0.859 


2.670 


-0.144 


Michaj 4 


34,403 


0.922 


2.081 


-0.004 


83,015 


0.947 


1.958 


-0.015 


42,145 


0.856 


2.698 


-0.102 


Michaj 5 


34,220 


0.922 


2.075 


-0.003 


82,273 


0.948 


1.949 


-0.017 


42,010 


0.857 


2.694 


-0.103 


Korf 1 


33,451 


0.924 


2.052 


0.115 


82,924 


0.947 


1.957 


-0.005 


41,766 


0.858 


2.686 


-0.063 


Korf 2 


33,620 


0.924 


2.057 


0.116 


82,328 


0.947 


1.950 


0.043 


42,758 


0.854 


2.717 


-0.081 


Korf 3 


33,638 


0.924 


2.058 


0.052 


81,940 


0.948 


1.945 


-0.017 


41,614 


0.858 


2.681 


-0.228 


Korf 4 


33,590 


0.924 


2.056 


0.117 


82,381 


0.947 


1.950 


0.051 


41,676 


0.858 


2.683 


-0.061 


Korf 5 


35,581 


0.924 


2.116 


0.117 


81,191 


0.948 


1.936 


0.057 


41,805 


0.858 


2.687 


-0.075 


Korf 6 


33,494 


0.924 


2.054 


0.109 


81,158 


0.948 


1.936 


0.063 


41,134 


0.860 


2.666 


-0.056 
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Figure 2. Residuals of Korf model 1 for spruce along tree 
diameter 
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Figure 3. Residuals of Korf model 1 for silver fir along tree 
diameter 
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Figure 4. Residuals of Korf model 1 for beech along tree 
diameter 

not for larger trees as expected. In addition, it is im- 
portant that there is no trend in residuals along DBH, 
which makes the function applicable. Although the 
bias was small for all variants, it is important that the 



bias is normally distributed throughout the diameter 
range. Therefore, mean residuals were calculated for 
the 10 cm diameter classes (tab. 3-5). The residuals 
were most appropriately distributed for Korf function 
(1) and Michailoff function (5). Silver fir (tab. 4) had 
the lowest residuals at ±0.1 m, most closely approxi- 
mated to the normal distribution throughout the entire 
diameter range. Higher bias was found only in trees 
exceeding 70 cm in diameter, and spruce (tab. 3) and 
beech (tab. 5) had slightly higher bias than fir trees. 
Meanwhile, Michailoff’s model (5) had similar, but 
somewhat negatively biased residuals for all species. 
This was most likely due to less flexibility in the mod- 
els based on Michailoff’s formula. 




Figure 5. Height-diameter models for spruce, fir, and beech 
based on Korf function 1 applied to LTP 45 on the first 
measurement when the stand was 82 years old, and on the 
last measurement when the stand age was 123 years, (dg and 
hg are mean stand quadratic diameter and height) 

Although all 11 models had relatively similar values 
for the estimated parameters, according to bias distribu- 
tion we propose that the optimal model is Korf 1. This 
was derived from Korf function (4) and contains 3 pa- 
rameters, A, k and n. Parameter A was fixed using d q , 
and h q , while k and n were dependent on the mean stand 
diameter and height d q , h q according to the product of 
their power (7.1) Since the height-diameter models were 
constructed for mixed-species forests of spruce, fir and 
beech with wide diameter and height variability, they 
were compared using the individual LTP’s, rather than 
simulation stands where all species would have the same 
mean stand heights and diameters. The height-based 
dominance of spruce at 82 years of age, and especially 
41 years later, is clearly evident for LTP 45 (fig. 5). Al- 
though the position and the shape of the models of two 
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growth stages were symmetrical for spruce and fir, the in mixed-species stands. Since the models had quite 
curves for beech intersected each other. This was most similar “goodness of fit” characteristics, it may eventu- 

likely due to beech height-growth dynamics, especially ate, that a more optimal model than Korf 1 will be found 



Table 3. Bias (m) of the height-diameter models for spruce at the individual diameter classes (gray cell indicates the best model 
- Korf 1) 



N 1 


DBH 

(cm) 


Model Korf 1-6 


Model Michailoff 1-5 


Smelko 

(1987) 


1 


2 


3 


4 


5 


6 


1 


2 


3 


4 


5 


172 


7.6 


-0.10 


-0.02 


-0.16 


0.05 


0.07 


0.14 


0.21 


0.54 


0.11 


0.37 


-0.01 


0.20 


744 


15.9 


0.06 


0.10 


-0.02 


-0.02 


0.00 


-0.05 


-0.44 


-0.51 


-0.31 


-0.39 


-0.33 


-0.48 


1,607 


25.4 


0.01 


0.02 


-0.04 


0.02 


0.04 


0.01 


-0.09 


-0.15 


-0.08 


-0.09 


-0.03 


-0.02 


2,246 


35.1 


0.20 


0.20 


0.15 


0.23 


0.23 


0.23 


0.26 


0.26 


0.20 


0.23 


0.22 


0.37 


1,751 


44.8 


0.15 


0.15 


0.08 


0.15 


0.14 


0.13 


0.10 


0.13 


0.07 


0.08 


0.05 


0.15 


957 


54.3 


0.04 


0.03 


-0.03 


0.01 


0.00 


0.00 


-0.18 


-0.15 


-0.17 


-0.18 


-0.20 


-0.24 


327 


64.0 


0.47 


0.44 


0.41 


0.43 


0.40 


0.45 


-0.02 


-0.02 


0.10 


0.03 


0.04 


-0.27 


102 


74.1 


-0.10 


-0.19 


-0.16 


-0.13 


-0.15 


-0.03 


-0.91 


-0.96 


-0.70 


-0.76 


-0.68 


-1.40 


37 


84.7 


-0.70 


-0.93 


-0.80 


-0.73 


-0.79 


-0.62 


-1.96 


-2.10 


-1.75 


-1.67 


-1.49 


-2.70 


5 


92.1 


-1.25 


-1.66 


-1.39 


-1.32 


-1.42 


-1.25 


-2.95 


-3.19 


-2.73 


-2.49 


-2.20 


-3.90 


2 


104.2 


1.10 


0.66 


1.06 


1.18 


1.10 


1.48 


-0.89 


-1.18 


-0.59 


-0.31 


0.07 


-2.14 


7,950 2 


37.2 


0.12 


0.12 


0.05 


0.12 


0.12 


0.11 


0.00 


-0.01 


-0.01 


0.00 


0.00 


0.02 



Note: 1 - number of height measurements; 2 - last row indicates the sums for columns and overall averages. 



Table 4. Bias (m) of the height-diameter models for fir at the individual diameter classes (gray cell indicates the best model 
- Korf 1) 



N 1 


DBH 

(cm) 


Model Korf 1-6 


Model Michailoff 1-5 


Smelko 

(1987) 


1 


2 


3 


4 


5 


6 


1 


2 


3 


4 


5 


1,799 


7.5 


-0.03 


0.06 


-0.01 


0.22 


0.02 


0.06 


0.02 


0.09 


-0.08 


0.18 


-0.02 


0.57 


4,705 


15.2 


-0.10 


-0.02 


-0.14 


-0.11 


-0.01 


-0.01 


-0.17 


-0.21 


-0.10 


-0.20 


-0.10 


-0.25 


5,382 


25.1 


0.03 


0.10 


0.05 


0.07 


0.14 


0.14 


0.08 


0.06 


0.07 


0.05 


0.06 


0.00 


4,584 


34.9 


0.05 


0.09 


0.04 


0.15 


0.13 


0.14 


0.11 


0.11 


0.03 


0.09 


0.06 


0.11 


2,873 


44.5 


0.02 


0.05 


-0.02 


0.15 


0.07 


0.08 


0.03 


0.05 


-0.06 


0.03 


-0.02 


0.08 


1,444 


54.4 


-0.06 


-0.07 


-0.11 


-0.03 


-0.10 


-0.10 


-0.17 


-0.15 


-0.19 


-0.15 


-0.18 


-0.13 


583 


64.2 


0.05 


-0.05 


0.03 


-0.12 


-0.12 


-0.13 


-0.22 


-0.23 


-0.11 


-0.18 


-0.14 


-0.25 


227 


74.0 


0.10 


-0.10 


0.12 


-0.32 


-0.19 


-0.22 


-0.37 


-0.43 


-0.11 


-0.29 


-0.15 


-0.50 


45 


84.1 


-0.45 


-0.70 


-0.35 


-1.09 


-0.87 


-0.89 


-1.08 


-1.10 


-0.63 


-0.94 


-0.80 


-1.30 


11 


93.9 


-0.64 


-0.94 


-0.65 


-1.51 


-1.24 


-1.25 


-1.46 


-1.46 


-1.06 


-1.28 


-1.11 


-1.75 


6 


105.7 


0.55 


0.00 


0.61 


-0.75 


-0.22 


-0.29 


-0.65 


-0.80 


0.07 


-0.41 


0.00 


-1.12 


2 


116.5 


1.12 


0.55 


0.84 


0.52 


1.01 


0.88 


-0.11 


-0.50 


0.08 


0.06 


0.65 


-0.45 


21,661 2 


29.8 


-0.01 


0.04 


-0.02 


0.05 


0.06 


0.06 


-0.01 


-0.02 


-0.03 


-0.02 


-0.02 


0.00 



Note: 1 - number of height measurements; 2 - last row indicates the sums for columns and overall averages. 
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Table 5. Bias (m) of the height-diameter models for beech at the individual diameter classes (gray cell indicates the best model 
- Korf 1) 



N 1 


DBH 

(cm) 


Model Korf 1-6 


Model Michailoff 1-5 


Smelko 

(1987) 


1 


2 


3 


4 


5 


6 


1 


2 


3 


4 


5 


446 


8.5 


0.50 


0.57 


0.52 


0.38 


0.30 


0.56 


0.43 


0.39 


0.71 


0.33 


0.38 


-1.13 


2,188 


15.0 


-0.34 


-0.31 


-0.54 


-0.32 


-0.31 


-0.36 


-0.38 


-0.36 


-0.47 


-0.35 


-0.36 


-1.57 


1,589 


24.7 


0.10 


0.06 


-0.12 


0.10 


0.09 


0.10 


0.07 


0.07 


-0.10 


0.08 


0.08 


-0.08 


878 


34.3 


0.16 


0.06 


-0.03 


0.12 


0.09 


0.15 


0.13 


0.11 


0.07 


0.10 


0.11 


0.78 


417 


44.3 


-0.12 


-0.20 


-0.24 


-0.13 


-0.16 


-0.07 


-0.17 


-0.19 


-0.07 


-0.21 


-0.19 


1.23 


142 


54.3 


-0.07 


-0.12 


-0.05 


0.00 


0.02 


0.00 


-0.18 


-0.14 


0.18 


-0.18 


-0.20 


1.82 


80 


63.8 


-0.46 


-0.57 


-0.49 


-0.34 


-0.33 


-0.38 


-0.56 


-0.53 


-0.13 


-0.56 


-0.57 


1.91 


36 


74.8 


-0.60 


-0.60 


-0.24 


-0.30 


-0.27 


-0.57 


-0.51 


-0.46 


0.12 


-0.42 


-0.37 


2.04 


12 


83.3 


-0.86 


-0.86 


-0.42 


-0.48 


-0.40 


-0.91 


-0.78 


-0.66 


0.09 


-0.61 


-0.65 


2.10 


6 


93.2 


-0.87 


-1.04 


-1.19 


-0.39 


-0.01 


-0.52 


-1.83 


-1.38 


-0.18 


-1.54 


-1.93 


3.18 


5,794 2 


24.4 


-0.06 


-0.08 


-0.23 


-0.06 


-0.08 


-0.06 


-0.10 


-0.10 


-0.14 


-0.10 


-0.10 


-0.40 



Note: 1 - number of height measurements; 2 - last row indicates the sums for columns and overall averages. 



using different datasets. However, if this occurred there 
is a high probability that such a model would be part 
of the Korf function-family. Although the models using 
Michailoff function were relatively satisfactory, they 
were less flexible for mixed-species stands with wide 
diameter and height variability, because they have sim- 
pler shape and less number of regression parameters to 
be estimated. 

Regression functions of the height-diameter curves 
express tree height depending on their diameter, mean 
stand diameter and mean stand height. Although the re- 
gression parameters were estimated with high precision 
and the explained variability was rather high (85-95%), 
there is still question if there are other factors (for in- 
stance social position of trees, stand density, stand mix- 
ture, etc.) that would significantly contribute to explain- 
ing the variability in tree heights. 

Comparison between the height-diameter 
models 

The models derived in this study were not evaluated us- 
ing an independent dataset. They were, however, com- 
pared to the existing models developed by Smelko et 
al. (1987) and derived using the graphical height-diam- 
eter curves published by Halaj (1955). For model com- 
parison and to avoid possible systematic errors in our 
model, Smelko’s model was tested by the empirical ma- 



terial used in this study. Residuals between measured 
height and those predicted by the Smelko model were 
calculated, together with bias for average tree diameter 
(tab. 3-5) and RMSE. Biases (m) for fir (tab. 4), spruce 
(tab. 3), and beech (tab. 5) were 0.004, 0.018 and -0.404 
and RMSE (m) were 1.99, 2.11 and 3.06, respectively. 
It is evident from the average bias that the models had 
systematic errors for all species. For spruce and fir, the 
predicted heights were systematically higher than the 
measured ones for tree diameters exceeding 55 cm. For 
beech, the predicted heights were higher when diam- 
eters were less than 25 cm and they were lower in larger 
sized trees, where differences ranged from 2-3 m in the 
largest trees. 

Temesgen and Gadow (2004) derived 10 height-di- 
ameter models for each of 8 species in British Columbia 
on permanent sample plots. The models were based on 
the three-parameter function (Yang et al. 1978). Two pa- 
rameters involved independent variables such as the ba- 
sal area of larger trees, the number of trees per hectare 
and basal area per hectare. The lowest RMSE (m) 1.44 
was found for cedar, but other species’ models recorded 
RMSE from 1.95 to 2.44. The spruce model had RMSE 
of 2.34 m and bias of 0.101 m, and therefore higher val- 
ues than our results. The original models developed for 
spruce in British Columbia were not subjected to testing 
for our dataset as it was done for the model derived by 
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Smelko et al. (1987) because it used different independ- 
ent variables for model generalization. 

Discussion 

It is generally known that the tree height depends not 
only on the diameter, age, species and site, but also on 
the mixture of species in the stand. Most of the stud- 
ies that have been conducted so far have proposed such 
models for poor forests, and most of growth models 
developed for Central European conditions included 
the relationships in poor forests: BWIN (Nagel 1999, 
2009), Moses (Hasenauer 1994; Kindermann and Ha- 
senauer 2005), Prognaus (Hasenauer and Monserud 
1996; Monserud and Sterba 1996; Nachtmann 2006), 
Silva (Pretzsch 1992; Kahn 1995) and SIBYLA (Fab- 
rika 2005). Only recently, research has been focusing 
on mixed-species forests, and several studies proposed 
height-diameter models for mixed-species forests (e.g. 
Temesgen and Gadow 2004). In our study, we devel- 
oped height-diameter models for spruce-fir-beech for- 
ests in the Western Carpathians. However, the three 
species did not occur on each plot, and the forests were 
mixed with different combination of the species. Co- 
occurrence of the all species was found only on a half of 
the plots. The co-occurrence of fir and spruce was found 
also on a half of the plots. Despite these shortcomings 
the dataset includes stands with the age from 30 to 200 
years. Another advantage is that trees were repeatedly 
measured within the 40-45-year period. This help in 
evaluating outliers or imprecise height measurements of 
individual trees (Curtis 1967; Pretzsch 2009). 

We used own experiences as well as published 
knowledge for the selection of the mathematical formu- 
las. From among many, only two were identified as most 
appropriate for the study: two-parameter Michajloff’s 
(1943) and three-parameter Korf’s (1939) function. The 
criterion that the curve passes through the one point de- 
fined by mean diameter and height was implemented for 
the study. The advantage of this approach compared to 
the approach where the point woidd be defined by the top 
diameter and height was a higher precision and a lower 
bias in tree volume calculation. After applying the crite- 
rion, the number of regression parameters was reduced 
to only one (Michajloff 1943) or two (Korf 1939). If the 
parameters are expressed depending on some stand 



characteristics, it defines a final shape and position of 
the model curve. Tree diameter, mean stand diameter 
and height are the parameters that rather precisely de- 
termine the tree height. Stand age as the other independ- 
ent parameter did not improved significantly the perfor- 
mance of the height-diameter models. An influence of 
other factors, for example top stand diameter and height 
(Hui and Gadow 1993; Stankova and Dieguez-Aranda 
2013) would probably not contribute to increase the pre- 
cision of the models. The influence of the stand density 
and stand mixture on tree height can be, however, dis- 
cussed. Temesgen and Gadow (2004) and Stankova and 
Dieguez-Aranda (2013) expressed the stand density by 
the number of trees per hectare. They based their models 
on the knowledge that the trees growing in more dense 
stands have usually higher heights. In our case, when the 
stands are consisted of both broadleaved and coniferous 
species, intra-species competition is compensated by 
inter-species competition. 

Conclusions 

Tree height is a very important variable characterizing 
stand structure. It provides an essential basis for tree 
volume calculation and site quality assessment. Al- 
though it has rather wide variability, it can be defined 
by a particular height-diameter model at a specific stand 
growth stage. Height-diameter models of one or several 
stands can be organized in a mathematical model which 
generalizes the growth of trees and forest stands at par- 
ticular stages throughout their lifespan. The general- 
ized height-diameter models were derived from a great 
amount of experimental data, including height measure- 
ments on 50 long-term experimental plots established 
in spruce-fir-beech mixed-species forest stands. Five 
models based on Michailoff’s function and six mod- 
els based on Korf’s function were tested for each spe- 
cies. Each model was adjusted so that the curves were 
constrained to pass through the mean stand height and 
diameter. Although tree height varied greatly within 
stands, all models had rather similar “goodness of fit” 
characteristics. The lowest RMSE errors of 1.94-1.97 m 
were found for silver fir, while spruce had 2.05-2.12 m 
and beech registered 2.67-2.72 m. These findings cor- 
respond to biological knowledge of the studied species. 
The final generalized model explained 95% of height 
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variability in silver fir, 92% in spruce, and 85-86% in 
beech. The remaining unexplained variability is pre- 
sumed to be due to factors not included in the models 
for various reasons. In these mixed-species forests, this 
may have been due to the horizontal and vertical distri- 
bution of the trees. However, it is difficult to explain this 
effect of horizontal and vertical arrangement in these 
three ecologically different species, not only at the time 
of measurement but also prior to measurement. 
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